LIFSCHITZ TAILS AND LOCALISATION 
FOR A CLASS OF SCHRODINGER OPERATORS 
WITH RANDOM BREATHER-TYPE POTENTIAL 



WERNER KIRSCH AND IVAN VESELIC 

f*"^ i Abstract. We derive bounds on the integrated density of states of Schrodinger operators 

with a random, ergodic potential. The potential depends on a sequence of random variables, 
not necessarily in a linear way. An example of such a random Schrodinger operator is 
the breather model, as introduced by Combes, Hislop and Mourre. For these models we 
■ show that the the integrated density of states near the bottom of the spectrum behaves 

X^/*X ' according to the so called Lifshitz asymptotics. This enables us to prove localisation in 

, certain energy /disorder regimes. 



1. Introduction 



> 

| In this paper we study spectral properties of certain Schrodinger operators with ran- 

■ dom potential. The spectral theory of such operators has been studied since the eighties 

in the mathematical literature and there are several monographs devoted to this topic, see 
e.g. [CFKS87, CL90, PF92, StoOl]. Certain spectral features, like the non-randomness of the 
i-C ■ spectral components and the integrated density of states, are shared by a wide variety of 

^p" 1 ! models under mild ergodiciy and regularity assumptions. However specific characteristics - 

43 ' like the existence of a certain spectral type — depend on the concrete model at hand. 

Our aim is to establish for a class random Schrodinger operators the Lifshitz asymptotics of 
the integrated density of states and localisation near the bottom of the spectrum. To explain 
this in more detail, we will define the considered class of operators and thereafter present our 
• i— i , results precisely. 

X 

Our results concern random Schrodinger operators H w of the following type. Let W per 
be a Z rf -periodic potential whose positive part V pei:+ := max(0, Vp er ) belongs to L\ QC and 
whose negative part V pcr - := max(0, — Vper) is in the Kato class. Let H w = Ho + W w and 
Hq = —A + W pcr , where is a random potential, i.e. a stochastic field, given by 

(1) W u {x):= J2< X k(u),x-k). 



Here A& : Q —> [\-, \+],k £ Z rf is a collection of non-trivial, independent, identically dis- 
tributed random variables on the probability space (fi, J-, P). The distribution of Ao is de- 
noted by \x and we assume inf supp fj, = A— . The function u is called single site potential and 
enjoys throughout the paper the following properties. 
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Assumption 1. The single site potential u: lxl d — > R is called monotone in the randomness 
if it is jointly measurable and satisfies the following: 
For all A G [A_, A + ] we have 

(2) supp^OcAx^h^]" 

d d 

(3) — u(A, •) G L°°(Ai) with k\ := sup sup —u(X,x) < oo 

"A xm d AG[A_,A+] OA 

For all x G R d and A G [A_, A + ] we have 

(4) a\ (A ' x) ^° 

There exist ei, 62 > such that for all A G [A_, A_ + £2] we have 

< 5 > i 



l dx u(X, x) G [ei, e 1 1 ] 



For all A G [A_, A+] the function 

(6) u(X, •) is reflection invariant with respect to all d coordinate axes. 

Note that the randomness enters the potential (1) via a field of random variables A&, k G Z d , 
not necessarily in a linear way. Condition (6) can be dispensed with. For this it is necessary to 
use in Section 3 instead of Neumann boundary conditions certain mixed boundary conditions 
introduced by Mezincescu in [Mez87]. Furthermore, one can assume instead of suppn(A, •) C 
Ai that there is a compact subset of M. d which conatins the support of u(X, •) for every A. A 

random potential of the form (1) with a single site potential satisfying Assumption 1 gives 
rise to a metrically transitive or ergodic operator, see e.g. [Kir89] or [PF92] for the definition. 
In particular, there is a subset £ of the real line such that the spectrum of H w concides with 
£ almost surely. 

Since A_ = infsupp^/ one can use Wely sequences to see that a(H per ) = Eq, where Eq 
denotes the the minimum of the almost sure spectrum of the family H u ,uj G il. 



Example 1. If we set in (1) u(A,x — k) = Xf(x — k) we obtain an alloy type potential 

(7) W u {x) := Y, f( x ~ k ) 

kez d 

Such random potentials have been thoroughly studied before in the context of the Lifshitz 
asymptotics of the integrated density of states and localisation, see e.g. [KM83, KS86, Mez87, 
CH94, Sto99, StoOl, GK01, GK04, AEN+06]. If / is non-negative and sufficiently regular the 
resulting single site potential u is monotone in the randomness. For such alloy type models 
the results we are aiming at are by now well understood, therefore we will not elaborate on 
them further. 



Example 2. More interesting is the situation, if we set 
(8) u(X,x) = -f(Xx) 



LIFSCHITZ TAILS FOR THE BREATHER MODEL 



3 



In this case the resulting stochastic field 

(9) WUx):=-J2 f^kH(x-k)) 

kez d 

is called random breather-type potential, cf. [CHM96, CHN01]. If we assume for the function 
/ 

(10) supp/cA A _, /€C o 1 (R d \{0}) 

(11) L°°(R d ) 3 g(x) := -x • (V/)(x) > and not identically vanishing 

(12) / is reflection symmetric with respect to all coordinate axes. 

then the potential u: R x M d — > R is monotone in the randomness. Inequality (11) is called 
the repulsivity property of /. 

Set 

(13) H pci := H + V pei and V pei (x) := ^ u(A_,x - k) 

kez d 

Remark 3 (Monotonicity of the model). Condition (4) ensures that for all oj G fl,x G M. d ,k G 
Z d we have 

(14) u(Xk,x) > u(\-,x) and consequently W w (x) > V pci (x) 

Denote for any A G [A_,A+] the modified single site potential u(X, x) := u(X,x) — u(\-,x) 
and 

V u (x) := W u (x) - V pei (x) = u(Xk(u),x- k) > 

hence H w = H pcT + V u . 

To formulate our first result we introduce some more notation. Let A L (j) := [-L/2,L/2] d + 
j C R d be a cube of side length L centered at j G Z d , II := Al n Z d the lattice points 
contained in it, and Xj the characteristic function of Ai(j'). We write A^ for A^(0). For a 
Schrodinger operator H on M d we denote by H i 'P er ) H L,N , H L,D its restriction to A^ with 
periodic, Neumann, respectively Dirichlet boundary conditions. Since we will be mostly using 
Neumann b.c. be abbreviate H L,N further to H L . Let xj(H^) denote the spectral projection 
for the operator H% associated with a Borel set Jet. 

An important spectral characteristic of ergodic random operators is the integrated density 
of states (IDS) denoted by N. It is defined as the limit of the distribution functions 

(15) N%(E) := L~ d n-th eigenvalue of H% does not exceed E} 

= L- d Tv[ X] _^ E] (H^)]. 

as L tends to infinity. Under our assumptions on H^^uj G £1, the limit exists at all continuity 
points of E and is independent of uj, for almost all wEd. If we replace in (15) the Neumann 
b.c. by periodic or Dirichlet ones, the limit distribution N does not change. Furthermore 
inf{£ | N(E) > 0} = E . 

Now we are able to formulate the result on the asymptotic behaviour of the IDS at the 
bottom of the spectrum. 
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Theorem 4 (Lifshitz Tails). Let H u ,uj £ Q be a random operator with potential (1) satisfying 
Assumption 1. Then 

(16) lim WWW < -t 

e\e log(E-Eo) - 2 



One can also prove equality in the above formula (16) if one assumes not to fast decay 
of fi([X-,X- + e]) as e — ► . We are only interested in the upper bound (16), which is the 
"hard" part of the equality, since we want to deduce spectral localisation. The bound (16) 
means that for E \ Eq asymptotically we have N(E) < ce~ cE d 

For breather type potentials satisfying Assumption 2 below we are furthermore able to 
prove localisation near the bottom of the spectrum. 

Let us explain this notion. An interval J = [Eq, Eq + 6] C R, 5 > is called a localisation 
interval for the family H^, u; £ 12, if has no continuous spectrum in J and all eigenfunctions 
associated to eigenvalues in J are decaying exponentially, for almost all u. In this situation one 
speaks of spectral or exponential localisation. Alternatively, localisation can also be expressed 
in terms of the dynamics of wavepackets, in which case one speaks of dynamical localisation. 
For a detailed discussion of this notion we refer to [StoOl] or [GK04]. Although dynamical and 
spectral localisation are not equivalent, for the type of models considered here it turns out that 
as soon as the multiscale analysis applies, both versions of localisation hold, see. e.g. [DS01, 
StoOl, GK01, GK04]. 

To carry through the multiscale proof of localisation one needs to establish beforehand 
several a priori conditions. The first one is the so-called initial length scale decay estimate on 
the integral kernel of the resolvent of a finite box Hamiltonian. It can be derived from the 
Lifshitz tail asymptotics of the IDS derived in our Theorem 4. This is a well know argument 
established e.g. in [MH84, Klo95, KSS98, StoOl, Ves02]. 

The second ingredient is a Wegner estimate, which goes back to upper bounds on the 
density of states established in the paper [Weg81]. Such estimates are well known for alloy 
type random Schrodinger operators as considered in Example 1. Some of the first papers 
where such estimates have been established are [Klo95, CH94, Kir96], see also [Ves04] for a 
survey. In [CHM96, CHN01] it was shown that Wegner estimates hold for random breather- 
type potentials, as long as the following assumption holds. 

Assumption 2. There exists an €4 > such that the single site potential /: R d — ► R of the 
breather-type potential satisfies for all x £ M. d 

(17) L°°(R d ) 3 -x ■ V/(x) > e 4 f(x) > 

The third assumption for the multiscale analysis independence of the potential values at a 
sufficiently large distance. This fact holds for the potentials considered in the present paper 
for the following reason: The single site potential has support in the cube Ai, for all values 
of A G [A_,A+]. Thus for two points x,y <G R d with \\x — y\\oo := maxf =1 \xi — yi\ > 1 the 
random variables uj 1— > V u (x) and uj \— ► V ul (y) are independent. 

If a Schrodinger operator with random ergodic potential satisfies the above three conditions, 
the multiscale analysis yields spectral and dynamical localisation in a (small) energy interval 
containing the bottom of the spectrum mfa(H u ), see e.g. [DS99, DBGOO, StoOl, GK01, 
GK04]. 
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Consequently one obtains the following result on localisation. 

Corollary 5 (Localisation). Let H W: uj G Q be a random operator with breather-type potential 
(1) satisfying Assumptions 1 and 2. Then there is a 5 > such that [E , Eq+S\ *s a localisation 
interval for H w ,lo G fl. 

Let us sketch the strategy of proof of Theorem 4 which is carried out in Section 2. 

To establish Lifshitz asymptotics of the integrated density of states we need to analyse 
the properties of eigenvalues close to the minimum of the spectrum and their associated 
eigenfunctions. More precisely, we 

(i) reduce the problem to bounds on the position of the lowest eigenvalue of the random 
Schrodinger operator. This is possible since we are interested only in the leading expo- 
nent in the relation (16). 

(ii) identify the configuration of random variables {\k\kel d such that the corresponding 
potential produces the minimal (with respect to the randomness) ground state energy 

(hi) establish how close a random potential has to be to the minimal potential configuration 
to produce an eigenvalue £ close to the overall minimum of the ground state energies. 

(iv) apply a large deviations estimate to bound the probability that such "close by" config- 
urations of the potential occur. 

Actually the above steps are implemented for restrictions of the considered Schrodinger op- 
erators on finite cubes with increasing length scales L ~ £~ 1 / 2 . Moreover on each scale we 
have to modify the random variable by a scale dependent cut-off. 

2. Proof of Lifshitz Tails 

Denote by ip\ the L 2 -normalised, positive ground state of Hp^ CT and by ^ its periodic 
extension on the whole of M. d . Then tpL '■= L~ d l 2 ~fy is the normalised ground state of Hp^ ev . 
Denote by E\(Hp^ cv ) the lowest eigenvalue of Hpef er . Using singular sequences and the 
fact that a pointwise positive L 2 -eigenfunction must be the ground state we conclude that 
mma(H pci ) = E^H^D = E L (H^f cr ) for all L G N. 

Remark 6 (Spectral gap for the periodic operator). Condition (6) ensures that ^> l is also 
the ground state of H P c T for all L G N. Thus min a(H pcr ) = Ei(Hp CT ) = E^Hper )• From 
[KS87] we infer that there exists a constant eo > such that 

E 2 (H^) - E 1 (H^) > £0 L- 2 

for all L G N. By adding a constant to the periodic potential we may assume without loss of 
generality that min a(Hp Cr ) = 0. Thus 

(18) E 2 {Hp£) > e L- 2 for all L G N 

Since the ground state \I> is pointwise uniformly bounded away from zero C3 := inf^gjjd ^(x) > 
0. We abbreviate by da(x) the measure ^(x) 2 dx on R d and use C4 := sup^Rd \I>(x) > 0. 

Since the dependence of the potential on the coupling constants Xk,k is not linear we 
introduce new random variables k. The scalar product we are interested in can be expressed 
as a sum in terms of these new random variables. 
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Remark 7 (Mapped random variables). Introduce for a parameter c 2 < f^-, i. e. < 
the cut-off random variables 

A fc := min{A fc , A_ +c 2 L -2 } € [A_, A_ + c 2 L~ 2 ] 

and the non-linearly mapped random variables 

£fcM = £(Afc) := y da(x)u(A fc (u;),x - k) 
Denote := H pcT + K, and K, := Y^k&A u(\k{oj),x - k). 



Remark 8 (Analysis of the first moment). The following expectation value will play a crucial 
role in the sequel 

(19) (Vl, hZ' N 1>l) = H^iPl) + (^l, Vu' n iI>l) = (V>l, K^l) 

For c 2 L~ 2 < €2, i. e. Ao < A_ + e 2 , we have 

£(Ao) = J da(x)u(Xo, x) — j da(x)u(\-,x) = J ^t— J da(x)u(r,x) 

< ^^c 4 by (5) 

ei 

< c 2l^. by definition of Ao 

ei L l 

HencG 



and thus for v := ^-L 2 + (ipL, H^' N ipL) we have 
by the choice of of C2 and (18). 

Remark 9 (Analysis of the second moment). We will need also an estimate for the second 
moment. By a zeroth order Taylor expansion one sees that for some A € [A_, A] 

~2 M a ,» \ wx ,du(X,x) 



u 2 (\,x) = 2(A- A_)u(A,x)- 



9A 

By (4) we have < u(A, x) < u(A, x) and thus 



("77 / /\ 

•u 2 (A, x) < 2(A - A_)u(A, x) — (A, x) 



< 2(A - A_)u(A, x) sup —(A, x) < 2k 1 (X - A_)«(A, x) 
xgR d oA 
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by Assumption (3). Hence 



/ 



da(x)u 2 (\k, x — k) < 2kic 2 L 2 j da(x)u(Xk, x — k) = 2k,\c 2 L 2 £fc 



and 

(21) \\Ht N ^L\\ 2 - L- d I da(x)u 2 (\ k ,x - k) < 2 Kl c 2 L' 2 L- d £ k 

kei L kei L 

The next theorem provides us with a lower bound on the first eigenvalue of a random box 
Hamiltonian. It is formulated in terms of an average of the random variables k. To prove it 
we use the Temple-inequality. The bounds on the first and second moment derived in the two 
preceding remarks are used on one hand to show that Temple's inequality is at all applicable, 
and on the other hand to insert them into Temple's inequality to obtain an appropriate lower 
bound. 

Theorem 10. Choose c 2 small enough such that 4k ^ 2 < \. Then 

kei L 

Proof. To ensure that Temple 's inequality can be applied to the operator H^ ,N and the vector 
V>L, we need to establish a chain of inequalities, see for instance Theorem XIII. 5 [RS78]. 
Assume c 2 L~ 2 < 62, then 

= Ek{H%£) < Ex(H^ N ) by monotonicity (4) 

< (tpL, H^' N ipi) by the min-max Theorem 

< v since eo > 

< E 2 (H^) by inequality (20) 

< E 2 (H^' N ) by monotonicity (4) 

We have checked the prerequisites for Temple's inequality and may apply it to the operator 
H^' N and the vector ipL: 

II U L ' N ,I, 1 1 2 

El {Ht N ) > {^H^l) ' ^ til , ^ 



2K lC2 L- 2 L- d Y, ke iM 



> ty L ,Hi? N ij, L ) 



to 



ui t-li/ eo t —2 

~2 L 



60 kTi L 

Here we used equation (21). 

It follows E 1 (H^ N ) > (l - ) (iPl,Hu' N iPl) = (l - ) Ekeijki") and thus we 
have proven the Theorem. □ 

The theorem im turn implies an estimate on how low most of the random variables k £ II 
are lying, if the principal eigenvalue of H% is low. 
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Corollary 11. Let < ±, c 6 > and set £ = j% > 0. ITten we /iave 

E 1 (H^ N )<£ implies #{£; G / L | £ fe < 45 } > ^ 
Proof. If the conclusion is false then 

€ II I 6 > 45} > ^ 

Hence Efce/ L >^ I T = 2 ^ Ld - lt follows from Theorem 10 

Ei(H^ N ) > \h- d ■ 28 L d = ^£> E 1 (Ht N ) 
which yields a contradiction. □ 
Now we have to show that the event 

#{k e I L | & < 45} > ^ 

has an exponentially small probability in the parameter L rf . To this aim we transform back 
first to the random variables A&, k and then to A^, k. 

Lemma 12. For c-j > 4/ei and C75 < £2, i.e. 5 small enough, we have 

£fc G [0,45[ implies G [A_,A_ + cy5[ 

Proof. Assume Afc > A_ + C75. By (5) we have for cj£ < 62 

> y da(x)u(X- + cj£,x) — J da(x)u(\-,x) 
r\-+c 7 £ d 



= / ^ T ~r~ / da(x)u(r, x) 

J\- dr J 



> c-j£ei 

> 45 since C7 > 4/ei 

It follows > 45 which is a contradiction. □ 

Now choose C6 < ^l. Then Afc G [A_, A_ + C7C6-L -2 ] implies A& G [A_, A_ + C7C6-L -2 ]. Thus, 
for c 7 > 4/ei, 2f < e ? 2 and c 6 < f^, 

E!(H^ N ) < g implies G J L | A fe < A_ + c 7 c 6 L- 2 } > ^ 

Now standard large deviations results imply 

Lemma 13. Let {^k}kez d be collection of independent, identically distributed random vari- 
ables on the probability space (f2,P) with E{Ao} > £3. Then there is a constant c§ such that 

¥{#{k G II I Xk < e 3 } >Y } ~ e ~ CSLd 
For suffciently large L we have E{Aq} > A_ + c-jc^L^ 2 and can thus apply the lemma. 
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